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The Gaussian transition in the spin-one Heisenberg chain with single-ion anisotropy is extremely 
difficult to treat, both analytically and numerically. We introduce an improved DMRG procedure 
with strict error control, which we use to access very large systems. By considering the bulk entropy, 
we determine the Gaussian transition point to 4-digit accuracy, Dc/J = 0.96845(8), resolving a 
long-standing debate in quantum magnetism. With this value, we obtain high-precision data for the 
critical behavior of quantities including the ground-state energy, gap, and transverse string-order 
parameter, and for the critical exponent, v = 1.472(2). Applying our improved technique at J- = 0.5 
highlights essential differences in critical behavior along the Gaussian transition line. 
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The Gaussian transition appears in several fields of 
quantum physics and statistical mechanics. The equiva- 
lence between surface-roughening transitions in classical 
two-dimensional (2D) models and quantum phase transi- 
tions in spin chains was introduced in Ref. [l|, and their 
rich phase diagrams investigated at length in Ref. y. 
Characterized by continuously variable exponents, the 
Gaussian transition differs significantly both from regular 
phase transitions and from those of Kosterlitz-Thouless 
(KT) type. These differences complicate both analytical 
and numerical approaches to a complete and accurate 
description of rough surfaces and quantum spin chains. 

The S = \ Heisenberg chain is one of the fundamental 
models in quantum magnetism. It formed the basis of 
Haldane's conjecture |3| for a finite gap in antiferromag- 
netic chains with integer spin, as opposed to the gapless 
spectrum of half-odd-integer cases. Numerically, quan- 
tum spin chains are important test-cases for any com- 
putational technique, and Haldane's prediction has been 
verified by a range of methods with increasing accuracy 
[J, |5[ . Experimentally, while the "Haldane gap" has been 
found in the excitation spectra of many systems [6[ , most 
known S = 1 chains, including NENP 0, NINAZ Q, 
and NDMAP [9|, are organic Ni materials with signifi- 
cant single-ion anisotropies. Analytical approaches to the 
Gaussian transition driven by this term are coniplicated 
by the lack of a suitable effective field theory [10[, and 
its broad nature makes all numerical techniques difficult 
to apply. Many authors have considered this transition, 
producing occasionally contradictory results lll-[l9j . 

In this Letter we resolve the problem of the Gaussian 
transition in the S = 1 chain with single-ion anisotropy. 
We exploit the fact that this transition is a gapless point 
between two gapped phases, whence the entropy exhibits 
a sharp peak. We introduce an improved density-matrix 
renormalization-group (DMRG) approach with system- 
atic error control, allowing high-precision calculations at 



system sizes up to L = 20000, which automatically elim- 
inate the end-spin entropy. We determine the critical 
point with very high accuracy, and thereby deduce the 
critical behavior of several quantities at different points 
on the Gaussian transition line. 
The general form of the model is 

^ = Eti^(^r^m+^f^^i) + J^sts!^, + D {stf (1) 

where J^ interpolates between XY and Ising spins, D is 
the single-ion anisotropy, and L the length of the chain. 
The full parameter space of {D,Jz) contains Neel, Hal- 
dane, large-!?, ferromagnetic, and two XY phases. In 
classical planar surface, or "solid-on-solid," models, the 
Neel and large- Z? phases are different "flat" phases, the 
Haldane phase is "rough," and the Gaussian transition 
is of "preroughening" type. These are the three phases 
of the S = 1 Heisenberg chain (J^ = 1) as D is var- 
ied. While the Neel phase possesses Z2 symmetry and 
the Haldane phase an incomplete Z2 x Z2 symmetry, the 
large- 13 phase has no remaining symmetries. The Gaus- 
sian transition is a line in the {D, J^) plane, on which the 
excitations are gapless. This line is well described by a 
conformal field theory (CFT) [20| , and has been analyzed 
in a number of studies 15l4l8l| , but none has achieved the 
numerical precision required for a consistent discussion of 
the critical behavior across the transition. 

DMRG is the most efficient and accurate numerical 
technique for ID systems [5|. Anticipating the need for 
both large system sizes and extreme precision, we be- 
gin by introducing an improved DMRG technique. In 
the conventional scheme, the absolute (coupled round-off 
and truncation) error increases systematically with L, 
and this accumulated error has a strong effect on the re- 
lability of the computation, possibly even disguising the 
critical behavior in a quantum many-body system. We fix 
the round-off error by renormalizing the lowest eigenvalue 
of Hamiltonian to remain of order 1, thereby obtaining a 



very significant reduction in the truncation error for large 
systems. In tlie DMRG iteration, we replace the original 
Hamiltonian matrix H(rn,L), for chains of L sites with 
m kept states, by H{m,L) — [ei{m,L — 2) — S], where 
El (m, L — 2) is the lowest eigenvalue of H{m, L — 2) and S 
is a constant chosen such that £1(771, L) ^ 0{1). Here we 
use S/J = 1 throughout. While the (extensive) total en- 
ergy of the ground state, Eg{m, L), can be reconstructed 
by summation, its (intensive) average value per site is 
determined directly and self-consistently as eg{m,L) = 
[ei(m, L) — S]/2. Similarly, for the first excited state 
ef{m, L) — [e2{m, L) — £2(771, L — 2)]/2 + eg{m, L), where 
62(171, L) is the second-lowest eigenvalue of H{m,L). 

The gap in our method is given simply by A(?Ti, L) = 
e2{m,L) — ei{m,L). Its general expression is 



A{m,L) = [ef{rn) - eg(rn)]L + A(rn) +X],'^i 






(2) 



where Cg (m) and e/ (m) are the intensive energies of the 
ground and first excited states for infinite L, and be- 
come equal for infinite m. In the polynomial expansion 
of contributions at higher order in 1/i, the 71 = 1 term 
arises from truncation errors and open boundary condi- 
tions (OBCs), while the 7i = 2 term has contributions 
from fluctuations at the quadratic band minimum. Here 
we calculate the energies in the linear term independently 
by extrapolation. Subtracting these gives a gap function 
A{m, L) that decreases monotonically with increasing L. 
A second polynomial fit of A (771) allows its extrapolation 
to infinite m to obtain the true gap. 

Sharing its foundations with quantum information the- 
ory, the DMRG method is ideally suited to discussions of 
entropy and entanglement. The von Neumann entropy, 
S{m,L) = — Trp(r7i,L) Inp(r77,-L), is readily computed 
from the reduced density matrix, which we obtain to high 
accuracy throughout our calculations with the renormal- 
ized Hamiltonian. The entropy obeys an area law except 
in critical regimes, where it depends logarithmically on L 
[21| . This extremum in entropy is an excellent indicator 
of a (gapless) critical point between two gapped phases. 

Before analyzing the entropy, we discuss the special 
and remarkable feature of the S = \ Heisenberg chain, 
that free S = 1/2 entities are found at a chain end, both 
in theory and in experiment [22|. In the Haldane phase 
with OBCs, the two free end-spins can be described by 
HeS = JbsSl'Sr [5[, where the effective coupling Jcs > 
falls exponentially with L. In the Hilbert space 5"*°' = 0, 
the two spins are maximally entangled with entropy In 2, 
while for 5'*°* = 1 they are unentanglcd. The additional 
truncation error due to this edge-entropy contribution 
causes significant computational difficulties. In Fig. 1(a) 
we find a In 2 drop in the ground-state entropy S(L) in 
the Hilbert space 5"*°' = when the chain reaches a cer- 
tain length at fixed D. For D = 0.92 and m = 1000, 
this occurs at L = 4500 [inset, Fig. 1(a)]. When L 
becomes sufficiently large, J^s falls below the machine 
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FIG. 1: (color online) Entropy <S as a function of D for J^ — 1. 
(a) Calculations with m = 1000. Open symbols are obtained 
for the lowest energy level in Hilbert space S'°' = with a 
range of L values, solid symbols for S*°* = 1. Inset: In 2 drop 
in S{L) for D = 0.92. (b) Bulk entropy 5(D) close to the 
Gaussian critical point, computed with L — 10000 for a range 
of m values. Insets: fitting slopes Al.i and Ar,i (left axis) 
and transition D^' (right) obtained as functions of m. 



precision and the end-spin contribution vanishes. The 
remaining "bulk" entropy contains the essential physics 
of the spin chain. Different but conceptually similar ap- 
proaches have considered both the two-site entropy and 
S{L) in a chain with no end-spin effects [23[. 

Figure 1(a) contrasts the total and the bulk entropy. 
Calculations with small L cannot access the unentangled 
regime, and for larger L we find a In 2 jump when D 
approaches Dc- For L = 10000. the end-spins remain 
entangled for 0.94 < D < Dc- The maximum in the 
total entropy moves str ong ly with L, showing no direct 
indication of criticality [17[. By contrast, in the Hilbert 



space iS*^" = 1, the end-spins are unentangled 



the 



lowest-energy state and this data reproduces exactly the 
bulk entropy. The location of the maximum in S, shown 
in detail in Fig. 1(b), is clearly invariant with m. A linear 
fit to the bulk entropy on both sides of the transition in 
Fig. 1(b) gives our primary result, Dc/J = 0.96845 with 
a minuscule error bar of 0.00008. The increasing slopes 
of the bulk entropy lines as 777 — > cx3 [inset Fig. 1(b)] 
indicate the onset of critical behavior. 

Having determined this extremely precise value of Dc, 
we may now discuss the critical behavior of the Gaussian 
transition with hitherto unattainable accuracy. We con- 
sider the physical quantities used in previous analyses of 
the transition 
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■19|, beginning with the gap. To avoid 
effects in the gap extrapolation related to the disappear- 
ance of edge states, we use the lowest energy levels in the 
Hilbert spaces 5'*°' = 1 and 5"*°' = 2. Figure 2(a) illus- 
trates our two-step extrapolation approach to compute 
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FIG. 2: (color online) (a) Gap as a function of L, computed 
for D I J = 0.95 with several values of m. Solid symbols for 
L -^ oo are extrapolated to m — > oo (inset), giving A(0.95) — 
0.00130(4). (b) Extrapolated gaps as a function of \D — Dc\- 



the gap for the extremely numerically challenging point 
D = 0.95, which lies very close to Dc- By following this 
procedure for all values of D, we show in Fig. 2(b) the 
approach of the gap to zero at Dc from both the Haldane 
and large-D sides. The closest four points, D = 0.925, 
0.95, 1.0, and 1.025, reveal a very narrow critical region, 
\D - Dc\ < 0.1, with critical exponent i/ = 1.472(4). 

In a CFT for the Gaussian critical line l20j , the gap A 
varies linearly and the energy eg quadratically with 1/L. 
For the CFT analysis, we perform DMRG calculations 
with periodic BCs (PBCs) using L = 200 and m = 2000 
[Figs. 3(a) and (b)]. We obtain the ground-state en- 
ergy eg = -0.86856650(4)J, velocity v = 2.564(2)J, 
central charge c = GfH/nv — 1.0006(8), Luttinger pa- 
rameter K ~ v/Aa — 1.321(1), and critical exponent 
V = 1/(2 -A') = 1.472(2). This last agrees exactly with 
our gap data in Fig. 2(b), confirming the consistency 
and accuracy of our calculations. Our computed cen- 
tral c harg e is extraordinarily close to the expected value 
c = 1 15|, ll6| . Even data at the extreme precision we at- 
tain cannot determine whether the second derivative of 



eg has a discontinuity [Fig. 3(c)], but set a very low upper 
bound. A continuous function with a point of inflection 
at Dc is consistent with the CFT expectation [l7| that 
the Gaussian transition be third-order for J^ = 1. 
The transverse string-order parameter is defined as 

O [l) ^ (5J exp (z^EU^p) Sf) , (3) 

and encapsulates the incomplete Z2 x Z2 symmetry of the 
Haldane phase [l[. To reduce the complexities inherent 
in calculating this quantity, we compute correlation func- 
tions only far from the system boundaries j24| , in the left- 
central block [i/4 - 1000, L/4] of the chain. We take the 



FIG. 3: (color online) (a) Finite-size extrapolation of low- 
est two gaps at Dc- Fitting lines from CFT give A{ST^ = 
0) = 2-Kv/L with V = 2.564(2)J and A(5'l;°* = 1) = 2-kJ(x/L 
with a = 0.48516(1). (b) Extrapolated ground-state energy 
eg at I?c, with CFT fit -0.86856650(4) J - /3J/L^ and /3 = 
1.3429(4). (c) Second derivative of extrapolated energy, (d) 
Extrapolated transverse string-order parameter (see text) of 
the Haldane phase, with fitting line 0.6036(4)|D - D^|0-353(i) 
Calculations for (a) and (b) performed with PBCs and m = 
2000, for (c) and (d) with OBCs and m = 1000. 



5*°* = 1 sector as the ground state. Figure 3(d) shows 
the results of our extrapolations to infinite L and m. 
The string-order parameter clearly shows excellent scal- 
ing behavior in the critical regime. The scaling exponent 
v' ~ 0.353(1) is very close to the value l/-\/8 predicted 
in the 2D classical model [1| , demonstrating the common 
physics of the Gaussian, or preroughening, transition. 

We illustrate with one example the utility of our 
improved DMRG calculations for investigating the en- 
tire Gaussian transition line. The point J^ = 0.5 has 



been considered by several authors [15l - ll8{ . Our results 
(Fig. 4) provide the most accurate information yet avail- 
able for this transition: Dd J ~ 0.6355(6). The values of 
L required to approach criticality are very much larger 
than for J^ = 1 [Fig. 4(a)], and the accuracy is lower be- 
cause S{D) is a significantly flatter function [Fig. 4(b)]. 
Our calculations with PBCs give eg = -0.91510889(1) J, 
V = 2.185(2)J, c = 1.000(1), K = 1.581(1), and v = 
2.387(5) at Dc, allowing a complete characterization of 
the physics of continuously varying exponents. 

We have considered the entropy 5(to, L) at finite m 
and L. In fact our results in Fig. 1 for tti = 1000 and L = 
10000 are fully converged for all values of D outside the 
very narrow region 0.94 < D < 1.00. We can deduce the 
critical behavior of S around Dc from a massive quantum 
field theory [U, in which S = (c/6)ln^ + A with ^ = 
v/A the correlation length and A oc \D — Dc^ ■ The 
convergent behavior of our data near Dc gives exactly 
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FIG. 4: (color online) S{D) as in Fig. 1 for J^ = 0.5. (a) 
Values of L as indicated. Inset: In 2 drop in S{L) for D = 0.4. 
(b) Values of m as indicated. 

the critical form S = Sq — {ciy / 6) In \D — Del, which is 
shown as the sohd lines diverging at Dc in Figs. 1 and 4. 

The Gaussian transition in the S ~ 1 chain is topolog- 
ical, in that the parity of the ground state changes from 
negative in the Haldane phase to positive in the large-D 
phase. The transition is thus associated with a change 
in the topological spin Berry phase from tt to [25| , and 
can be followed by a method of crossing energy levels 
(of states in the appropriate parity sectors). Our high- 
precision results demonstrate that this is indeed a very 
sensitive indicator of a topological transition: among all 
previous studies llMl9l|. we find that the only accurate 
estimate of Dc was obtained, despite being limited to 
16-site systems, by employing this approach lla |. 

We have demonstrated that the entropy is very valu- 
able for discussing continuous phase transitions between 
gapped states. Many other types of strongly interacting 
quantum system fall in this category, one good exam- 
ple with electronic de gree s of freedom being the ionic 
Hubbard model (IHM) [26[ . The numerically challenging 
transition in this case is of KT type. Continuous gapped- 
to-gapped transitions for both bosonic and fermionic sys- 
tems exist in ultra-cold atomic condensates on optical 
lattices. The Gaussian transition has not yet been ob- 
served in experiment, due to difficulties in controlling the 
ratio D/J in condensed matter systems, and cold-atom 
experiments may offer a clean solution to this problem. 

To summarize, we calculate the critical point of the 
spin-one Heisenberg chain with single-ion anisotropy, 
Dc/J = 0.96845(8), to extremely high accuracy. To 
achieve this we introduce an improved DMRG scheme, 
which controls the absolute error of a large system and 
allows the elimination of end-spin effects. We exploit this 
accuracy to deduce the critical properties of many quanti- 
ties at the Gaussian transition. The energy, entropy, and 



gap all show good scaling behavior with a single critical 
exponent, ly = 1.472(2). We apply our technique also at 
Jz = 0.5 to illustrate the continuous variability of expo- 
nents on the Gaussian transition line. 
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